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ABSTRACT 

To maximize the number of planet detections, current microlensing follow-up observations are focusing on 
high-magnification events which have a higher chance of being perturbed by central caustics. In this paper, 
we investigate the properties of central caustics and the perturbations induced by them. We derive analytic 
expressions of the location, size, and shape of the central caustic as a function of the star-planet separation, s, 
and the planet/star mass ratio, q, under the planetary perturbative approximation and compare the results with 
those based on numerical computations. While it has been known that the size of the planetary caustic is oc yfq, 
we find from this work that the dependence of the size of the central caustic on q is linear, i.e., oc q, implying 
that the central caustic shrinks much more rapidly with the decrease of q compared to the planetary caustic. 
The central-caustic size depends also on the star-planet separation. If the size of the caustic is defined as the 
separation between the two cusps on the star-planet axis (horizontal width), we find that the dependence of the 
central-caustic size on the separation is oc (s+s -1 ). While the size of the central caustic depends both on s and q, 
its shape defined as the vertical/horizontal width ratio, 7Z C , is solely dependent on the planetary separation and 
we derive an analytic relation between 1Z C and s. Due to the smaller size of the central caustic combined with 
much more rapid decrease of its size with the decrease of q, the effect of finite source size on the perturbation 
induced by the central caustic is much more severe than the effect on the perturbation induced by the planetary 
caustic. As a result, we find that although giant planets with q > 10~ 3 can be detected from the planet search 
strategy of monitoring high-magnification events, detecting signals of Earth-mass planets with q ~ 10~ 5 will 
be very difficult. Although the central caustics of a pair of planets with separations s and s~ l are identical 
up to the linear order, we find that the magnification patterns induced by the pair of the degenerate caustics 
of planets with q > 10~ 3 are different to the level of being noticed from observations with < 2% photometry. 
Considering that the majority of planets to be detected by the strategy of monitoring high-magnification events 
are giant planets, we predict that the s «-* s~ l degeneracy could be broken for a majority of planetary events 
from observations with good enough precision. 

Subject headings: planetary systems - planets and satellites: general - gravitatinal lensing 



1. INTRODUCTION 

Various methods have been proposed to detect and 
characterize ext rasolar planets, i ncludi n g radial veloc- 
ity t echnique (iMavor & Ouelozl 119951 iMarcv & Butlerl 
1996), transit method jSjruyd 119521), direct imaging 
(Angej 11994 IStahl & Sandler! Il995l) . pulsar timing 
analysis ( Wolszczan & Fraill 119921). and microlensing 
jMao & Paczvnskil 1199 1[ [Gould & Loebl 119921) . See also 
the review of|Perrvman (2000). The microlensing signal 
of a planetary companion to Galactic disk and bulge mi- 
crolens stars is the short-duration perturbation to the smooth 
standard light curve of the primary-induced lensing event 
occurred on background star. Compared to other methods, 
the decay of the planetary lensing signal with the decrease 
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of the planet/star mass ratio is relatively slow, and thus the 
microlensing technique has an important advantage of being 
applicable to the detections of Ea rth-mass planets b y using 
already existing instrument (Ben nett & Rhid fl996). The 
microlensing method also has a uniqu e applicability to the 
detections of free-floating plan ets ( iBennett & Rhid 120021 
lHan et al.lf2004 lHan. et al.ll2005l) . In addition, the method is 
not restricted to planets o f nearby stars and can b e extended 
even to nearby galaxies (Baltz & Gond ololl2001l) . Recently, 
two clear- cut microlensing d etec tions of exoplanets were 
reported bv lBond et ail J2004I) and lUdalski et alJ (120051) . 

Due to the rare and incidental chance of lensing events 
combined with the short duration of planet signals, plane- 
tary lensing searches require a special observational setup, 
where survey observations issue alerts of ongoing events 
in the early stage of le nsing magnification (Soszvhs ki et alJ 
1200 It [Bond et all 1200 ll) and foll ow-up collaborat ions inten- 
sively monitor the a lerted events (B ond et alJ2002llPark et ail 
2004; iCassan et al] l2004). However, follow-up is generally 
done with small field-of-view instrument, and thus events are 
monitored sequentially. As a result, only a handful num- 
ber of events can be followed at any given time, limiting 
the number of planet detections. To maximize the num- 
ber of potential planet detections wi t h a lim ited use of re- 
sources and time, Griest & Safizadeh ( 1998) proposed to fo- 
cus on high-magnification events. They pointed out that high- 
magnification events have a dramatically higher chance of 
being perturbed by planets due to the existence of a 'cen- 
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teal' caustic (see § 2 for more detail). By adopting this 
proposal, current foll ow-up experiments are giving hig h pri- 
ority to these events lAlbrow et alJl2001t lR attenburvt l2002l 
iBond et alJl2002l lAbe et al J 120041 Uiang et al] 120041) . How- 
ever, little has been studied about the characteristics of cen- 
tral caustics and the perturbations induced by them (central 
perturbations). 

In this paper, we investigate the properties of central caus- 
tics and the perturbations induced by them. The layout of 
the paper is as follows. In § 2, we describe basic physics of 
planetary lensing. In § 3, we analytically investigate how the 
location, size, and shape of the central caustic vary depending 
on the mass ratio, q, and separation, s, of planets under the 
planetary perturbative approximation and compare the results 
with those based on numerical computations. In § 4, we then 
systematically inspect the patterns of central perturbations by 
constructing maps of magnification excess for planets with 
various values of s and q. From the constructed maps, we 
also examine the effect of finite source size and investigate 
the possible types of planets detectable for given source stars 
with various sizes. We summarize the results and conclude in 
§5. 

2. BASICS OF PLANETARY LENSING 

The lensing mapping from the lens plane to the source plane 
of -point masses with no external shear or convergence is 
described by the lens equation of 



nij/M 

z-Zlj' 



(1) 



where £ = £ + irj, zlj = Xlj + iyij, and z = x + iy are the com- 
plex notations of the source, lens, and image positions, re- 
spectively, z denotes the complex conjugate of z, ntj are the 
masses of the individual lens components, and M = £\ trij is 
the total mass of the lens system dWiti 19901) . Here all angles 
are normalized to the Einstein ring radius 9e of the total mass 
of the system, i.e, 



?e = 



4GM ( 1 1 



1 V2 



(2) 



where Z\ and £>s are the distances to the lens and source, re- 
spectively. For a single lens (N = 1), there exist two images 
with locations at u; ± = 0.5[u ± (w 2 +4) 1 / 2 ] and magnifications 
of A± = 0.5(A± 1), where u = \(-zl\ is the separation be- 
tween the lens and source. Then, the total magnification cor- 
responds to the sum of the magnifications of the individual 
images and it is related to u by 



A=A + +A_ 



(m 2 + 2) 
u(u 2 + 4) 1 / 2 ' 



(3) 



A planetary lensing is described by the formalism of a bi- 
nary (N = 2) lens with a very low-mass companion. For a 
binary lens, there are three or five images and the number of 
images changes by two as the source crosses a caustic. The 
caustics are important features of binary lensing and they rep- 
resent the set of source positions at which the magnification of 
a point source becomes infinite. The caustics of binary lens- 
ing form a single or multiple closed figures where each figure 
is composed of concave curves (fold caustics) that meet at 
cusps. For a planetary case, there exist two sets of discon- 
nected caustics: one 'central caustic' located close to the host 
star and one or two 'planetary caustics' depending on whether 



X : primary star 

+ : effective primary star position 
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FlG. 1. — Geometry of a central caustic. 



the planet lies outside or inside the Einstein ring. Since the 
central caustic lies close to the host star, the perturbation in- 
duced by the central caustic occurs close to the peak of lensing 
light curves of high-magnification events. For a star-planet 
system, it is known that the lensing behavior can be described 
by the perturbation approach due to the sma ll planet/star mass 
ratio llBozzall 19991 lAsadall200l 1X3120051) . In this case, the 
lensing equation is expressed as 
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where the position of the star is chosen as the coordinate ori- 
gin and z p represents the location of the planet. 

3. PROPERTIES OF CENTRAL CAUSTICS 

Under the planetary perturbative approximation (q <C 1 and 
| \z p | - 1 1 ^> q), the location of the central caustic can be ex- 
pressed in an analytic form of 
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), where the polar coordinates are centered at the 
position of the primary star and tj> represents the polar angle 
(see Figure[Q. With the notation for the star-planet separation 
of s= \z P \, the caustic position is expressed in an explicit form 
of 



s + s + 2(cos </)— 2cosc/j) 
(s + s _1 -2cosc/)) 



-q- 



2sin J 



(s + s 1 -2cos</j) ; 



(6) 



(7) 



The central caustic has an elongated asteroid shape with four 
cusps, where two of them are located on the star-planet axis 
and the other two are located off the axis. Cusps of a central 
caustic occur when d( c /d(f> = (i.e., when c/> = and ir for the 
on-axis cusps and when <fi = <fi c and 2ir - <fi c for the off-axis 
cusps), and thus their locations are found to be 
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FIG. 2. — Variations of the shape and size of the central caustic produced by planets with various separations from and mass ratios relative to the primary star. 
The left panels show the variations of the shape of the central caustic depending on the separation, while the right panels show the variation depending on the 
mass ratio. In all cases, coordinates are centered at the effective position of the host star and the planets are located on the left side. Note that the axis scales differ 
from one panel to another for better comparison of the caustic shapes. 




-3.0 -3.5 -4.0 -4.5 -5.0 

log q 

FIG. 3. — Size of the central caustic, as measured by the horizontal width, 
as a function of s and q. 
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FIG. 4. — Dependence of the central-caustic shape, as measured by the 
ratio between horizontal and vertical widths, on the planet-star separation. 
The computations are based on planets with q = 10 -4 . 



where 
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Then, the horizontal and vertical widths of the central caustic 
defined as the separations between the on- and off-axis cusps 
(see Figure^ are respectively 
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In Figure|2] we present central caustics produced by planets 
with various separations and mass ratios. From the figure and 
equations (|5} - d!2i . we find the following characteristics of 
the central caustic. 

1 . The size of the central caustic depends on both the sep- 
aration and mass ratio. The dependence of the central- 
caustic size on the mass ratio is linear, i.e., cx q, as 
shown in equations (|5} - Q and demonstrated in the 
right panels of Figure[2] By comparison, the size of the 
planetary caustic is proportional to ^fq. This implies 
that the central caustic shrinks much more rapidly with 
the decrease of the planet mass compared to the plane- 
tary caustic. If the caustic size is defined as the horizon- 
tal width, the size depends on the star-planet separation 
by A£ c cx (s + s -1 ) -2 . Then, the caustic size becomes 
maximum when s ~ 1 and decreases with the increase 
of \s — 1 1. In the limiting cases of a very wide-separation 
planet {s ^> 1) and a close-in planet (s <C 1), the depen- 
dencies are respectively 



fors > 1, 
fors< 1. 



(13) 



In Figure[3] we present the size of the central caustic as 
a function of s and q. 

2. For a given mass ratio, a pair of central caustics with 
separations s and s~ l are identical up to linear order. 
This can be seen from the pair of caustics with sepa- 
rations s and s" 1 presented in Fig. [2] Analytically, this 
can also be proved from equations (|3 - Q where the 



inversion of s - 



result in the same expression. 



3. Under the perturbative approximation, the shape of the 
central caustic is solely dependent on the planet sep- 
aration. We quantify the caustic shape as the verti- 
cal/horizontal width ratio, 1Z C = A?y c /A£ c , andpresent 
the variation of 1Z C as a function of s in Figure^ In the 
figure, we present two curves, where one is based on 
the perturbative approximation and the other is based 
on numerical computation. One finds that the width ra- 
tios based on the analytic and numerical computations 
match very well except the region around s ~ 1 , where 
the perturbative approximation is not valid. We find that 
the difference in this region is caused by the merge of 
the planetary and central caustics. 10 One also finds that 
the width ratio of the central caustic rapidly increases 
with the decrease of |j — 1|. 

4. CENTRAL CAUSTIC PERTURBATIONS 

Knowing now the characteristic of central caustics, we then 
investigate the pattern of planetary perturbations induced by 
central caustics. For this purpose, we construct maps of mag- 

10 The planetary caustic produced by a planet with a separation s from its 
parent star is located at x pc ~ s— 1/s. Therefore, as s — > 1, x pc — * 0, and thus 
the planetary caustic approaches the central caustic and eventually merge. 
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FIG. 5. — Contour maps of magnification excess as a function of source position (£,f?) for planetary lens systems with various star-planet separations and 
planet/star mass ratios in the region around central caustics. Contours are drawn at the levels of t = ±2% (thin curves) and ±5% (thick curves) and greyscale is 
used to represent positive (bright) and negative (dark) deviation regions. All lengths are scaled by the Einstein ring radius corresponding to the total mass of the 
lens system. The coordinates are centered at the effective position of the host star. In all cases, planets are located on the left side. The maps are constructed for 
a main-sequence source star with a normalized radius of p* = 0.0018, that corresponds to the case where the source with an absolute radius of R+ = 1 Rq and 
located at Ds = 8 kpc is lensed by a lens with m = 0.3 Mq and Dl = 6 kpc. 

FIG. 6. — Contour maps of magnification excess for the same planetary lens systems as in Fig. [3] but for events associated with a turn-off source star with 
p,= 0.0054. 



TABLE 1 

Detectability of Planetary Signals 



source 




detectability 




type 


q = 10~ 3 


q = 10~ 4 t 


j = 10~ 5 


main-sequence (p, = 0.0018) 


O 


O 


X 


turn-off (p* = 0.0054) 


O 


A 


X 


clump-giant (p* = 0.0234) 


o 


X 


X 



NOTE. — The adopted source radius are R+ = 1 Rq, 3 Rq, 
and 13 Rq for the main-sequence, turn-off, and clump giant stars, 
respectively. The normalized source radius p* = #*/<9e of each 
source star is determined by choosing an Einstein ring radius 
of £?e = 0.32 mas that corresponds to the value of the most fre- 
quent Galactic bulge event with m = 0.3 Mq, Z>l = 6 kpc, and 
Ds = 8 kpc. The individual symbols for the detectability represent 
possible ('O'), marginal ('A'), and very difficult ('X'), respec- 
tively. 



nification excess, which is defined as 
(t ^ A ~ A ° 



(14) 



where A is the exact magnification of the planetary lensing 
and Aq is the single lensing magnification caused by the host 
star at its 'effective' position. Due to the additional deflec- 
tion of light produced by the presence of the companion, it 
was known that the effective lensing position of a compo- 
nent of the binary lens system is shifted tow ard its companion 
JDi Stefano & MadH996l lAn & Hanll200l . The amount of 
the shift of a lens component T toward the other component 
7 is 

nij/nii 



I— >} 



(s + s- l )/(8 E j/9 E ) 8 E 



(15) 



where to, and mj are the masses of the individual lens com- 
po nents and Or, ■ i is th e Einstein ring radius corresponding to 
m, (I An & Hanl l2002). For a planetary lens case, the shift of 
the host star (with a notation '★') toward the planet (with a 
notation '/?') is expressed as 



q 



(s + s- 1 ) 1 



(16) 



because f?E,* ~ 8 E and trip/m* = q. u Then, as either |j— 1| — > 
oo (wide-separation planet) or \s— 1 1 — > (close-in planet), the 
shift is Axl^^j, — > 0, and thus the effective lensing position 
of the primary star approaches its original position. 

1 ' We note that the term ' i n equations I15t and II hi was not included 
in the corresponding equations of Di Stefano & Mao 1 1996) and An & Han 
1 2002), where they treated only wide-separation binaries. We introduce this 
term to keep the symmetry between the central caustics of planets with sepa- 
rations .v and y _1 . 



In Figure El we present the constructed magnification- 
excess maps for planets with various values of s and q. Plane- 
tary signal is important for planets with separa tions with in the 
so called 'lensing zone ' of 1/1.6 < s < 1.6 (Gould &Loebl 
cm, and thus we plot maps for planets located within this 
range. To see the effect of finite-source size, we construct 
three sets of maps with source stars of main-sequence (with 
a radius of R+ = 1 Rq), turn-off (with = 3 Rq), and clump 
giant (with R+ = 13 Rq) stars and they are presented in Fig- 
ure 13 |6] and respectively. The magnification of a finite 
source with a surface brightness profile of I(Q is computed by 
the intensity-weighted magnification averaged over the source 
star surface, i.e., 



Af s (0- 



J s I(C')A(C + C')dC' 



(17) 



where A denotes the point-source magnification, £ is the vec- 
tor position of the center of the source, is the displacement 
vector of a point on the source star surface with respect to the 
source star's center, and the two-dimensional integral is over 
the source-star surface S. The effect of finite source size is 
smearing out the detailed structures of the planetary lensing 
signals. The finite-source effect is determined by the normal- 
ized source radius p+, which represents the angular radius of a 
star, 0*, in units of the Einstein ring radius, i.e. p+ = 9+/8 E . We 
thus set of the individual stars by choosing an Einstein ra- 
dius of 8 E = 0.32 mas that corresponds to the value of the most 
frequent Galactic bulge event with to = 0.3 Mq, D e = 6 kpc 
and Ds = 8 kpc. Then, the normalized source radius of a 
R+ = 1 Rq star corresponds to = 0.0018. For a source with 
uniform surface brightness, the computation can be reduced 
from a two-dimensional to a one-d imensional integral using 
the Generalized Stokes's theorem (Gould & Gaucherel 1997; 
Dominik 1998). To accelerate computations, we thus assume 
that the source stars have uniform surface brightness. How- 
ever, we note that the effect of non-uniform surface brightness 
on the planetary lensing signal is not important. For a plane- 
tary system with a primary star of to = 0.3 Mq, the mass ratios 
of the Jupiter-, Saturn-, Neptune-, and Earth-mass planets cor- 
respond to q ~ 3 x 10~ 3 , 10~ 3 , 10~ 4 , and 10~ 5 , respectively. 

From the excess maps, we find the following properties of 
central perturbations. 

1 . First, we find that the effect of finite source size on cen- 
tral perturbations is much more severe than the effect 
on perturbations caused by planetary caustics. This is 
because the central caustic is not only much smaller 
but also shrinks much more rapidly with the decrease 
of q than the planetary caustic. As a result, we find 
that although giant planets with q > 10~ 3 can be de- 
tected from the planet search strategy of monitoring 
high-magnification events, detecting signals of Earth- 
mass planets with q ~ 10~ 5 would be very difficult. It is 



Fig. 7.— 
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Contour maps of magnification excess for the same planetary lens systems as in Fig.|5] but for events associated with a clump-giant with p* 



5 

= 0.0234. 



estimated that the lower mass limit of detectable plan- 
ets will be that of a Neptune-mass planet with q ~ 10 -4 , 
but detecting this mass-range planets will be possible 
only for events associated with source stars smaller than 
turn-off stars. In Tabled we summarize the possible 
types of planets detectable from events involved with 
various types of source stars. 

2. Second, as expected by the close similarity between the 
pair of caustics with separations s and s~ l , the pattern 
of perturbations induced by these two caustics are sim- 
ilar each other. A good example is the pair of maps 
with (s,q) = (1.2, 1(T 4 ) and (1/1.2, 10" 4 ) in Figure |5] 
However, we find that as the planet mass increases, the 
difference between the perturbation patterns of the pair 
of planets increases. We find that the difference is large 
enough to be noticed for planets with q > 10~ 3 , c.f., the 
pair of maps of the planets with (s,q) = (1.2, 10~ 3 ) and 
(1/1.2, 10~ 3 ). Considering that the majority of plan- 
ets to be detected by the strategy of monitoring high- 
magnification events are giant planets, we predict that 
the s «-> l/s degeneracy could be broken for a major- 
ity of planetary events from observations with a < 2% 
photometric precision. 

5. SUMMARY AND CONCLUSION 

We investigated the properties of central caustics and the 
perturbations induced by them. Under the planetary pertur- 
bative approximation, we derived analytic expressions of the 
location, size, and shape of the central caustic as a function 
of the planet separation and mass ratio and compared the re- 



sults with those based on numerical computations. While it 
has been known that the size of the planetary caustic is pro- 
portional to yfq, we found from this work that the dependence 
of the central-caustic size on the mass ratio is linear. As a re- 
sult, the central caustic shrinks much more rapidly with the 
decrease of the planet mass compared to the planetary caus- 
tic. Due to the large finite-source effect caused by the smaller 
size of the central caustic combined with much more rapid 
decrease of its size with the decrease of q, we predict that 
although giant planets with q > 10~ 3 can be detected from 
the planet search strategy of monitoring high-magnification 
events, detecting signals of Earth-mass planets with q ~ 10~ 5 
would be very difficult. Although the central caustics of a pair 
of planets with separations s and s~ l are identical up to the 
linear order, we found that the difference between the mag- 
nification patterns induced by the pair of degenerate caustics 
of planets with q > 10~ 3 can be noticed at the level of ~ 2%. 
Because the majority of planets expected to be detected by 
the strategy of monitoring high-magnification events are gi- 



ant planets, we predict that the s 



degeneracy could be 



broken for a majority of planetary events from observations 
with good enough photometric precision. 
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